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1. Introduction 



The equation of state (EOS) of quantum chromodynamics (QCD) has been perused for a long 
time, not only because it is theoretically interesting but also for its growing importance in the rela- 
tivistic heavy-ion collision experiments. Almost two decades back, an operator method formalism 
was devised [jl|] to study the EOS of QCD numerically by using lattice gauge theory. Later it was 
found that computations on coarse lattices, using this method, gave negative pressure near the phase 
transition region. It was argued then, that this problem arises due the use of perturbative formulae 
for the various derivatives of the coupling. To cure this problem a new elegant method, the integral 
method, was proposed in [Q]. This method assumes the system to be homogeneous, bypassing the 
use of those derivatives. It uses another non-perturbatively determined derivative of the coupling, 
namely the QCD j3 -function. Unfortunately, the assumption of homogeneity does not hold at a 
first order phase transition. Moreover, the method also has other limitations like the pressure be- 
ing ill-defined below an arbitrarily chosen cut-off temperature and the fluctuation measures being 
obtained only by taking numerical derivatives, resulting in large errors. 

In view of this situation we propose ^ a new variant of the operator method which gives 
non-negative, well-defined pressure for the entire temperature range. This method can also be 
extended to the calculation of fluctuation measures [Qj, starting from the first principle. We choose 
the temporal lattice spacing to set the scale of the theory, in contrast to the choice of the spatial 
lattice spacing in the approach of [Q]. Thus, our method [^] could be called the t-favoured operator 
method and the method of ^ may be called the s-favoured operator method. 



2. Formalism 



As in most finite temperature lattice gauge theory formulation, our t-favoured formalism is 
defined on a 3 + 1 dimensional hypercubic, asymmetric lattice having different lattice spacings 
a z and a s in the temporal and spatial directions respectively. The temperature(r) and the spatial 
volume(V) are defined as T = (a x N z )~ l and V = (a s N s ) , N z and N s being the number of lattice 
sites in the temporal and spatial directions. In our t-favoured scheme we introduce the anisotropy 
parameter E, and the scale a by the relations — 



and 



a = a T 



(2.1) 



The partial derivatives with respect to temperature and volume can easily be written in terms of 
these variables — 



*4- 

dT 



v 5 H 



da 



and 



3 dE, 



(2.2) 



In finite temperature lattice gauge theory simulations the scale is set by the temperature T = 
1 /a t N t . In that sense our choice of scale a = a z corresponds to the situations in actual simulations, 
giving rise to the most natural procedure for the scale setting. On the other hand, in the case of the 
s-favoured scheme for any £ 7^ 1, the scale is set by a = a s and only in the £ — > 1 limit this 
natural choice of the scale emerges. 
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For a pure gauge SU (N c ) theory on an asymmetric lattice, the Wilson action and the corre- 
sponding partition function are defined as — 



S[U]=K s P s + K t P T , 



and 



3T(y,T) = J We-W, 



(2.3) 



where periodic boundary conditions are assumed. P s denotes the sum of spatial plaquettes over all 
lattice sites, and P t is the corresponding sum of mixed space-time plaquettes. The couplings in the 
temporal and spatial directions are given by K s = 2N C / £,g 2 , and K x = 2N C E, jg\ respectively. 

In the weak coupling limit, gj 2 's (i = s, x) can be expanded [Q] around their symmetric lattice 
value g~ 2 (a), 

gr\a^) = g -\a)+ Ci ^) + 0[g\a)), (2.4) 

with the condition c,-(£ = 1) = 0. The value of g 2 (a) depends on the identification of the scale a. 

We are now in a position to calculate different thermodynamic quantities. Let us first look 
at the energy density (e) and the pressure (P). Starting from the standard definitions of £ and P 



in terms of In 3f, operating the lattice derivatives of eq. ( [L2| ) on the partition function given in 



eq. ( |2.3| ) and making use of the relation mentioned in eq. ( ]2.4| ) one can derive the expressions for e 
and P. In the E, — ► 1 limit these expressions reduces to — 



D s ~Dr , , , N 

-4-^- (c' s D s + c' x D % ) 
8 

-J-^-(c' s D s + c' T D % ) 



B(a s ) 
2%a 2 



D s +D x 



and 



(2.5) 



Here the primes denote derivatives with respect to B(a s ) is the QCD j3 -function, with the usual 
definition of a s = g 1 {a)/An. An ultraviolet divergence in the values of the plaquettes is removed 
by a subtraction of the corresponding vacuum (T = 0) values to yield D ( - = (Pi) — (Pq) above. In 
[Bj, c\ 's have been calculated for a weakly coupled SU(N C ) gauge theory upto one-loop order. 
Since these Karsch coefficients are known only upto one-loop order, we also use one-loop order 
B(a s ) for logical consistency. We use the one-loop order renormalized coupling g 2 (a) as has been 
suggested in [^]. 



Our expressions for energy density and pressure, eq. (|2.5D, can be compared to the expression 
of £ and P in the s-favoured operator method of [[I]]. It is clear that the P in the t-favoured method is 
exactly e/3 of the s-favoured method. The positivity of the £ in the s-favoured scheme guarantees 
the positivity of P in the t-favoured scheme. The interaction measure — 



A 

j4 



(g-3P) 

j4 



T 2na 2 



D s + D r 



(2.6) 



is same for both the t and s-favoured methods. Since A is always positive our expression for e is 
also bound to give positive values for the energy density. 

Let us now turn our attention to specific heat at constant volume (c v ) and the isentropic speed 
of sound (c s ). These quantities are defined as — 



dT 



and 



c: 



dp 


dP 




' de 


de 









-1 



s/T 3 
c v /T 3 



(2.7) 
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Here we have used the definition of entropy density, s = (dS/dV)r = (e +P)/T, and the thermo- 
dynamic identity (dP/dT)y = (dS/dV)r, S being the total entropy. In [Q] it has been argued that 
the specific heat can most easily be obtained by working in terms of the dimensionless variable c €, 
the so-called conformal measure, — 



A 



and 



r = r 



(2.8) 



By doing some straight forward algebraic manipulations it is very easy to see that- 



T 3 



\P/T 4 J 



s r e 

f3 + 3 fl 



and 



P/T A 
7/f 1 



1 + 



Ts/T 4 



(2.9) 



In [J3j] , an expression for Y has been calculated in terms of the plaquettes and the co- variances of the 
plaquettes. It has also been shown in [Q] that c v /T 3 and c 2 s reaches their respective ideal gas values 
in the g — > limit. The expression for T contains the second derivatives of the Karsch coefficients 
C; 's. Numerical values of these second derivatives have already been evaluated in for an SU (3) 
gauge theory. Using all these we can determine the specific heat and the speed of sound directly 
from numerical simulations. 



3. Results 

We have performed simulations in the temperature range 0.9 <T /T c < 3, with temporal lattice 
sizes N t = S, 10, 12. Spatial lattice sizes were chosen to be 2N T + 2 for T < 2T C , and 3N t + 2 for 
3T C . For the zero temperature simulations, a minimum lattice size of 22 were chosen and then 
lattice sizes were scaled up with changes in lattice spacings. Typical number of sweeps used for 
the measurements of the plaquettes were around few hundred thousand in finite temperature runs, 
and around one hundred thousand in zero temperature simulations. The details of our simulations 
can be found in [J3j] . 

From the computations for the above mentioned three temporal lattices we performed contin- 
uum extrapolations by linear fits in c k 1 /N% at all temperatures. In Figure (|l]a) and Figure ([l]b) 
we compare the continuum limits of pressure and energy density for the t-favoured method [Q], 
s-favoured method ^ and the integral method 0]. The most important thing to be noted is that, 
unlike the old s-favoured operator method, the new t-favoured operator method gives positive pres- 
sure for all the temperatures, and also for all lattice spacings. The t-favoured method pressure 
agrees well with the integral method pressure in the high temperature (T > 2T C ) region and both 
differ from the ideal gas value by about 20%. In contrast, even in the continuum limit, and on A^ T 
as large as 8 — 12, the s-favoured method does not yield positive pressure near T c . On the other 
hand in the phase transition region the t-favoured method pressure shows a steeper rise than that of 
the integral method pressure. Compared to the integral method the energy density in the t-favoured 
method is harder near T c and agrees with that of the integral method for T > 2T C . This indicates a 
difference in the latent heat determined by the two methods. 

Part of the disagreement between the t-favoured and the integral method results can be traced 
back to the fact the integral method pressure is simply set to zero at T c and below this cut-off 
temperature the integral method pressure is ill-defined. Also the assumption of a homogeneous 
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1 1.5 2 2.5 3 3.5 1 1.5 2 2.5 3 3.5 

T/T c T/T c 



Figure 1: We show comparisons between the continuum results of different thermodynamic quantities for 
our new method (boxes), the old operator method (triangles) and the integral method (line). In panel (c) 
we show a comparison between our continuum results for c v / T 3 (open boxes) and continuum 4e / T 4 (filled 
boxes). The data for the integral method has been taken form Ref. [||. 



system breaks down at a first order phase transition that takes place for a pure SU (3) gauge theory. 
Along with these, the integral method results of [Q] shown here are for coarser lattices, uses non- 
perturbative /3 -function and different scheme to determine the renormalized coupling. We plan to 
investigate these issues in future, in view of all the differences in these two computations. 

We present the continuum results for the specific heat and the speed of sound in Figure (|I|c) and 
Figure (jljd) respectively. For T > 2T C , c v /T 3 is far from its ideal gas value, but is quite consistent 
with the prediction in conformal theories that c v /T 3 = 4e/T 4 . In the neighbourhood of T c , c v 
shows a peak-like structure. Whereas c s is consistent with its ideal gas value (in accordance with 
the prediction for a conformal theory) for T > 2T C , it decreases dramatically near T c . 

We find that the pressure in the t-favoured method agrees with that obtained from a dimension- 
ally reduced theory, matched with the 4-d theory perturbatively upto order g 6 ln(l/g) [|[|, almost 
all the way down to T c . This has been shown in Figure (||a). 

In Figure (||b) we compare our result for the entropy density s with that of the strongly coupled 
jV = 4 Super Yang-Mills (SYM) theory obtained by using the AdS/CFT correspondence [^]. We 
find that the deviation of s/sq, so being the entropy density for an ideal gluon gas, from 3/4 is in 
good agreement with the predicted results for a strongly coupled J/ = 4 SYM from T > l.5T c . 
This, along with our results for c v /T 3 and c 2 s , gives a hint that QCD behaves like a quasi-conformal 
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Figure 2: In panel (a) we compare the pressure obtained by t-favoured method (boxes), integral method 
(dotted line) and the g 6 ln(l /g) order perturbative expansion (solid line). The data for the integral method 
and the perturbative expansion are taken from Ref. [|[| and Ref. |S] respectively. The values of the T /A^ 
in Ref. [||] has been converted to T /T c using the T c /Ajjg quoted in TO. In panel (b) we show the deviation of 
s/sq from 3/4 (boxes) as a function of the t'Hooft coupling. The line shows the same quantity as predicted 
by the AdS/CFT correspondence [||]. 



theory for t' Hooft couplings g 2 N c < 9. 
4. Summary 

We presented a new variant of the operator method, namely the t-favoured scheme, for the 
calculation of the EOS of QCD. This methods gives positive pressure for all temperatures and all 
lattice spacings, even when the old operator method gives negative pressure. Using this new 
method we have computed the continuum limit of the energy density and pressure. We have also 
extended this method to determine the continuum limit of specific heat at constant volume and 
speed of sound, both above and below the transition temperature. To the best of our knowledge 
this is the first continuum results for these two quantities. We have also compared our results with 
the predictions of perturbation theory and also with the results for strongly coupled jV = 4 SYM. 
We find good agreement with both in certain ranges of temperature. 
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